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ABSTRACT 

We calculate the impact of the intergalactic medium (IGM) on the observed Lya 
lines emitted by galaxies in an ionised IGM at z ^ 4. Our model accounts for gas 
clumping in the IGM and for the fact that high-redshift galaxies reside in overdense 
regions, which causes the velocity field of the IGM to depart from the Hubble flow. 
The observed shape of the Lya line varies widely, with dependence on the intrinsic 
width and systemic velocity of the line, a galaxies star formation rate and the local 
extra-galactic UV-background. For large star formation rates and levels of the UV- 
background, absorption in the IGM does not result in a Lya line that is asymmetric as 
is common among known high-redshift Lya emitters. For models in which the lines do 
show the observed strong asymmetry, the IGM typically transmits only 10 — 30% of the 
Lya flux. The increase in the ionising background that accompanied the completion 
of reionisation barely increased the IGM transmission, which suggests that LAEs of 
comparable luminosity should not appear to be significantly dimmer prior to overlap. 
In this light, we briefly discuss the potential of Lya emitters as a probe into the epoch 
of reionisation. 

Key words: cosmology: theory-galaxies: high redshift-line: profiles-radiative trans- 
fer 



1 INTRODUCTION 

The Lya emission line is a p owerful tracer of you ng galaxies 
JPartridge fc Peeblesl 1967) and quasars (e.g. iFan et all 
120061 ) in the high redshift Universe. Sin ce the discovery 
of th e first high-redshift Lya emitters |Hu fc McMahonl 
1 19961 : ISteidel et all 1 19961 : iTrager et all Il997l . hereafter. 
LAEs), existing Lya surveys have detected sev eral hun- 
dred g alaxies out t o redsh i ft z = 7.0 (e.g. iHu et alj 
19981: iRhoads et al.1 120001: [Malhotra fc Rhoadsl 12002 



For example, the shape of the Lya emission line as- 
sociated with outflows produced by superwinds is typi- 
cally peaked on a frequency redshifted by u w relative to 
the true line center, where f w is the speed of the outflow- 
ing gas. Furthermore, the spectrum is typically asym met- 
ric and contains more photons redward of this peak (|Ahn 
2004 lHansen fc Ohll2006l : lTasitsiomi|[2006l ; IVerhamme et al. 
2006). On the other hand, the Lya emission line from 
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Hu et al.1 12003: iKodaira et al.1 120031 ; iDawson et al.1 1 200 
Stanwav et al. | l2004l ; iTran et alj|2004 iTaniguchi et al.ll2005l: 
Westra et al.l I2OO6I; iTapken et all 120061: Kashikawa , et al 



200e : IShimasaku et al.ll2006l : live et al.ll2006l : IStanwav et a! 
20071 ). and the number of these known LAEs is steadily 
growing. The total flux in the Lya line emerging from 
galaxies provides a measure of its star formation rate, while 
its exact shape contain information on the sources powering 
the emission and on the gas distribution and kinematics 
through which the Lya photons propagate. 
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galaxies surrounded by neutral collapsing gas clouds is ex- 
pected to be blueshifted relative to the true line center, 
by an amount that depends on the infall velocity and 
the total column density of neutral hydrogen in the in- 
falling gas ilZheng fc Miralda-Escudll2002l: lTasitsiomill2006l : 
iDiikstra et al.ll2006al : IVerhamme et al.ll2006l )~ 

It is well known that the total flux in the Lya line and 
its shape are affected by subsequent scattering in the IGM. 
To first order, the IGM at high redshift erases the blue half, 
and transmits the red half of the line, reducing the total Lya 
flux by a factor of 2. This is due to the fact that photons 
blueward of the Lya transition will redshift into the Lya 
resonance, where there is a large probability for the photon 
to be scattered out of the line of sight. Photons redward 
of the Lya transition will redshift further and further away 
from resonance, which reduces the probability that these 
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photons are scattered out of the line of sight. This simple 
picture would apply if the gas surrounding the galaxy is 
expanding with the Hubble flow. However, infall of IGM 
gas onto galaxies will occur out to distances well beyon d 
the virial radius (lLoeb fc Eisensteinl [19951 ; iBarkanal l2004h . 
Furthermore, this infalling gas is expected to be significantly 
overdense, with values as large as A ~ 20 times the mean 
density (just outside the virial radius §0- 

Infall may have strong implications for the im p act o f 
the IGM on the Lya emitted by galaxies. ISantosI (|2004l ) 
has calculated the impact of an infalling neutral (i.e. pre- 
reionisation) IGM on a galaxies Lya spectrum and found 
that the IGM typically erases a large fraction of the red side 
of the Lya line. Therefore, much more than half of the Lya 
is lost, with the IGM transmitting only T a ~ 10 -3 — 0.3 
of all Lya photons. This low IGM transmission is partly 
due to the IGM damping wing, and partly due to resonant 
absorption by neutral hydr o gen atoms in the i nfall region 
(see Fig. 8 of ISantosI liooi ). IfBarkana fc Loebl |2003h have 
found evidence for gas infall in the Lya spectra of one bright 
quasar at z = 4.8. Despite its large ionising flux, denser 
clumps in the infalling gas absorb observable amounts of 
Lya, red ucing the total Lya tr ansmission well below 1. In 
addition, Dij kstra et alj (|2006bh have shown that infall m ay 
also explain the o bserved spectrum of [steidel et a l. (2000)'s 
Lya 'Blob' LABI (|Wilman et al.ll2005l ) and th e peculiar line 
shape of a z = 5.8 Lya emitter observed by iBunker et al.l 
(2003). These observations suggest that if infall is common 
around known LAEs, then their intrinsic Lya luminosities, 
and consequently, their star-formation rates, may have been 
(significantly) underestimated in many cases, since the Lya 
transmission is generally not corrected for infall. 

In this paper we address the question: How does the 
IGM affect the shape and transmission of th e Lya emission 
line? Our approach is similar to the work of ISantosI (|2004l ) 
and our results are consistent whenever there is overlap. 
However, we expand upon this work in several important 
ways: (1.) we incorporate gas clumping in the IGM; (2.) we 
investigate the impact of clustering of galaxies on the IGM 
transmission; (3) we calculate the skewness parameter, Sw , 
which was recently introduced by Kashika wa et alj {2006). 
However, the biggest difference between this paper and pre- 
vious work is our focus on the reionised IGM: although it is 
generally understood that a neutra l IGM can strongly sup- 
press Lya emission from g alaxies (|Gunn fc Peterson! 1 19651 : 
lHaiman|[2002l : ISantosll2004l ), the impact of an ionised IGM 
has been virtually ignored. As discussed in § !5.5l our analysis 
is also relevant in scenarios in which large disjoint HII bub- 
bles exist prior to the end of reionisation dFurlanetto et all 
l2004al : I Wyithe fc Loebl l2004al ) . Note that we focus on the 
Lya line emitted by galaxies. For a brief discussion on the 
impact of the IGM on Lya lines emitted by quasars, the 
reader is referred to lDiikstra fc Wyithe I l|2006r ). The outline 
of this paper is as follows: We review the mechanism that 
creates Lya emission in star forming galaxies, and discuss 
its expected properties in §[2] We describe our model for the 
IGM in § In § [3] we present calculations of the modified 
shape of the Lya line. We discuss various implications of 
our results in § O before presenting conclusions in § [6] The 
parameters for the background cosmology used throughout 
this paper are Q. m = 0-24, Oy = 76, Q b = 0.044, h = 0.73 
and a 8 = 0.74 jSpergel et al.ll2006l) . 



2 LYa EMISSION FROM STAR FORMING 
GALAXIES 

2.1 The Lya Luminosity vs. Star Formation Rate 

In this section we review the relation between the star for- 
mation rate in galaxies and Lya luminosity, and discuss its 
uncertainties. The Lya flux from star forming galaxies orig- 
inates in the dense nebulae from which the stars form. Ion- 
ising photons produced by O stars, which are absorbed in 
the ne bulae, are conver ted into Lya in about 2 out of 3 
cases (jOsterbrockl 1 19891 . for case-B recombination, p 367). 
Under the assumption that all Lya photons escape from the 
galaxy into the IGM (sometimes referred to as the 'dust- 
uncorrected' Lya luminosity) , the total Lya luminosity of a 
galaxy can then be estimated from 

L Lya = 0.68/w a [l - /obc]Qh, (1) 

where hv a = 10.2 eV is the energy of a Lya photon, / csc 
is the escape fraction of ionising photons, and Qh is total 
luminosity of ionising photons. As Qh depends on the num- 
ber of O-stars, its value depends on the assumed IMF and 
metallicity of the gas from which the stars form. Under the 
assumption o f constan t star formation rate (as opposed to 
a star burst) . ISchaererl (|2003T ) has calculated Qh for several 
different IMFs, and for a range of metallicities Z (expressed 
in solar units Zq). For a Salpeter IMF with lower and upper 
mass limits of Mi — Mq and M u = IOOMq, respectively 

log 10 Qh = 53.8 + log 10 M, - 0.0029(9 + log 10 Z) 2 ' 5 , (2) 

in wh ich M* is the star formation rate in Mq /yr. ISchaererl 
(|2003h also calculates the time-evolution of Qh for an in- 
stantaneous star burst. In this case, Qh can be very high, 
reaching values of log 10 Qh (photons s~ Mq 1 ) ~ 47 — 48, 
during the first few Myr after the star burst. For the re- 
minder of this paper, we will assume the star formation 
does not take place in such a burst and we will use Eq {1} 
and Eq ^ to calculate I/Lya as a function of M t . Following 
Santos! (|2004l ). our models assume / csc = 0.1 ijlnoue et al.l 
20061 ) and Z = 0.05Zq, unless stated otherwise^ This low 
metallicity reflects that galaxies at z >4 — 5 live in Universe 
that is <jl Gyr old, and in which there has been less time 
for their gas to become enriched by metals and dust. 

2.2 The Shape+Width of the Intrinsic Lya Line 

For simplicity, the shape of the intrinsic Lya line is usually 
assumed to be Gaussian. For gas that is optically thin to Lya 
photons, a reasonable choice for the standar d deviation o f 
this Gaussian emission line is a a ~ 1.5u c h-c (|Santosll2003 ). 
However, to take this value for all galaxies assumes that 
each star forming, Lya emitting galaxy is a disk seen edge- 
on, which is only true for a small subset of LAEs. For disk 

1 For comparison, given solar metallicities (Z = Zq) the total 
output of ionising photons is related to the star form ation rate 
through log 10 Q H = 53.0 + log 10 M, [|Kennicuttlll998ll . This re- 
lation is often used in the literature and is a lower by factor of 

2 than the relation used in this paper. The uncertainties in the 
relation between I/Lyc* an d M t introduced by the gas metallicity 
and by the IMF are ignored for the reminder of this paper. Here, 
our focus is on the uncertainties introduced by the IGM. 
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galaxies seen face-on, the Lya line is narrower. To account 
for this inclination effect, the fiducial intrinsic Lya emission 
line should be lowered to ^ ~ 0.7 x l.Sticirc ~ Vcirc (where 
(| cosi|) ~ 0.7, and i is inclination angle of the galaxy). 

As already mentioned, the previous discussion applies 
when the gas inside the galaxy is optically thin to Lya pho- 
tons. Since the line-center Lya optical depth exceeds unity 
for hydrogen columns Nh <T0 13 cm -2 , the optically thin ap- 
proximation is most likely incorrect. Indeed, as mentioned in 
§ 15. li the bulk of Lya photons are believed to be emitted in- 
side dense nebulae in which star formation is occurring. For 
ionising photons to be converted into Lya, these nebulae 
must be optically thick to ionising radiation. This requires 
hydrogen to have a column of Nh J>10 17 cm -2 , making the 
line center Lya optical depth to J>10 4 . Resonant scattering 
of Lya photons through optically thick neutral hydrogen gas 
clouds may boost the width of the Lya emission line. How- 
ever, resonant scattering is not expected to boost the Lya 
line width significantly. 

For example, if we have Lya photons emitted in the 
center of a static uniform slab/sphere with line-center opti- 
cal depth ro, then resonant scattering causes the emerging 
spectrum to comprise two peaks, each separated from the 
line center by ~ 180(iV H /10 20 cm- 2 ) 1/3 ( T gaB /10 4 K) 1/6 km 
s" 1 , where Tgas is the gas temper ature l|Harringtonl 1 19731 ; 
lNeufeldlll990l ; lDiikstra et al.ll2006ah . The FWHM of the res- 
onant profile is ~ 2.5 times as large as this separation. 

To quantify the impact of resonant scattering on the 
width of the line, we convolve a Gaussian emission line 
with the reson ant profile derived by iHarringtonl (|l973l ) and 
iNeufeldl (|l990T ). For this example, the Gaussian has a stan- 
dard deviation of w C j rc = 160 km s . The associated phys- 
ical picture is that Lya photons emerge from the nebula 
in which they were created with a double peaked emission 
spectrum, after which they can freely escape from the galaxy. 
For simplicity, we assume that the total column of hydrogen, 
Nh, is the same for each nebula. In Figure [1] we show Lya 
emission lines for various values of Nh- The figure shows 
that broadening of the line is not important for Nh <J10 20 
cm~ 2 . The main reason is that in this example, the FWHM 
of the Gaussian is comparable to the FWHM of the reso- 
nant profile for Nh = 10 20 cm~ 2 , which is ~ 450 km s _1 . 
For larger columns the Lya line broadens and displays the 
double-peak of the resonant profile. 

The column densities of atomic hydrogen observed in 
local galax ies typically lie in the rang e N H = 10 20 - 10 21 
cm -2 (e.g. lBosm^ll98ll ; lKu"nth et alj|l998h . but may reach 
larger values. This suggests that resonant scattering can 
boost the observed Lya widths significantly. However, 
the above calculation represents an upper limit on the 
amount of resonant broadening for a given Nh, because the 
calculation assumed the hydrogen gas to be uniform and 
static. In reality, velocity gradients in the gas will facilitate 
the escape of Lya from the nebula, resulting in a reduced 
broadening through resonant scattering. Furthermore, the 
observation that ionising radiation does emerge from star 
forming regions suggests that paths through the galaxy 
must exist along which Nh <J10 17 cm -2 . If Lya photons 
escape through these lower column-density paths, then 
resonant broadening of the Lya line is negligible. This also 
applies when the ISM consists of two phases in which case 
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Figure 1. Broadening of the Lya emission line through resonant 
scattering. The Lya flux density J(x) (arbitrary units) is shown as 
a function of velocity off-set from the true line center. A Gaussian 
Lya emission line with u v = 160 km s~ x , was convolved with the 
'resonant profile', which is the line profile that would emerge from 
a static, uniform slab with a hydrogen column density TVjj (see 
text). When Nh j>10 20 cm -2 , resonant scattering of the Lya 
emission line becomes more important. 



resonant scattering is even less important (see § 12. 3f) . The 
possibility that the majority of Lya photons escape through 
these 'low column density paths', has been invoked as an 
explanation that t he Lya photons manage to esca pe from 
galaxies at all (e. dNeuferdlll99ll : iKunth et al.lll998T l. 

We conclude that resonant scattering of Lya photons is 
not expected to broaden the width of the Lya emission line 
of galaxies with u C i rc <;150 — 200 km s~ by a factor Jj2. 
For significantly less massive galaxies, i.e. v c i IC <J50 km s _1 , 
resonant scattering could in principle broaden the line to 
a few hundred km s _1 (since this requires columns of only 
Nh ~ 10 20 cm -2 ) though it may not if the Lya escapes 
through 'low column density paths'. 



2.3 The Equivalent Width of the Lya Line 

If we approximate the spectrum of a source blueward of the 
Lya frequency using a powerlaw of the form f v = ku 13 , then 
the m aximum intrinsic Lya EW is given bv (ICharlot fc Falll 
1 19931 . for (3 < -1) 



EW„ 



0.68his a J t 



oa f(y) 



dv 



A(A = A„) 



-0.68- 



3 UJ 



(3) 



where h = 6.62 x 10~ erg s, is Planck's constant and 
vh = 3.28 x 10 15 Hz, is the hydrogen ionisation threshold. 
The Lya flux was calculated assuming that two-thirds of all 
ionising photons are converted into Lya, i.e. we assumed the 
escape fraction of ionising photons to be / csc = 0. EW max 
is then obtained after dividing this number by the contin- 
uum flux density (in erg s _1 A -1 ) at the Lya frequency. 
Here, we used the relation A/a = vj v and have denoted the 
Lyman limit frequency with vh- For stars, the index is 
typically —4 to —5 and we get that EW max = 64 A (for 
(3 = -4). ICharlot fc Falll (119931 ) showed, using stellar syn- 
thesis models in combination with realistic stellar spectra, 
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that the maximum EW that is reached in the first ~ 10 7 
yrs after the onset of star formation and can be as high as 
EW m ax ~ 300 A, a lthough the pre cise value depends on the 
assumed IMF fe.g. lSchaererll2003l ). 

Resonant scattering of Lya photons increases the to- 
tal distance traveled before the photons emerge from the 
galaxy. Therefore, Lya is more subject to dust attenuation 
than the continuum photons, which can result in a reduc- 
tion of the Lya EW . However, the opposite may be tru e in 
a multi-phase ISM (|Neufeldl [l99ll ; H arisen fc Ohll2006l ). in 
which case the dust is locked up in cold clumps embedded 
in a hot ionised medium. The Lya photons can escape from 
the galaxy scattering off the dust clumps, while mainly trav- 
eling through the ionised interclump medium. Depending on 
the details of clump size, clump dust content, etc, the frac- 
tion of Lya transmitted through the ISM could exceed that 
of continuum photons. Therefore, the intrinsic Lya EW pro- 
duced by a normal stellar population could, in some cases, 
be larger than 300 A. 



MODELING THE INTERGALACTIC 
MEDIUM 



Our a pproach to modeling the IGM is based on iBarkanal 
who calculated the average density and velocity pro- 
files of intergalactic gas around virialised halos of various 
masses. It is worth emphasising that infall of gas onto galax- 
ies is inherent in the Press-Schechter model for the collapse 
of cold dark matter halos. This model is described below. 



3.1 The IGM Density and Velocity Field 

Figures 2 and 4 of Barkana (2004) show the average radial 
density and velocity profile of the intergalactic gas beyond 
the virial radius r v j r of objects of mass Mtot at various red- 
shifts. Here, we focus on the density profile associated with 
2 = 5 and M tot = 10 10 M Q - 1O 12 M . Just outside the virial 
radius, the average density, p(r) is p(r v i r ) = 20p, where p is 
the mean baryon density of the universe. The density drops 
off approximately following a power-law p(r) oc r _1 , and 
reaches a value of p(r = 10r v j r ) = 2p. At larger radii, the 
density asymptotically approaches p (see Fig 2 of Barkana, 
2004). The velocity profile (Fig 4 of Barkana, 2004), shows 
the IGM to be collapsing at ~ v c [ IC , just outside the virial 
radius. The IGM is stationary with respect to the galaxy 
at r ~ 5r v i r - In this model, the IGM never becomes fully 
comoving, which is not surprising since the average density 
within a sphere of radius r is always larger than p. The model 
should therefore break down beyond a certain radius. In this 
paper, we assume that this model applies up to 10r v j r . Be- 
yond this radius, we assume the IGM is simply comoving 
and at the mean density in the universe. To summarise, our 
fiducial IGM density and velocity fields are given by 



p(r) 



wigm(j") 



20p{r/r vir y 
P 



\10r vlI H(z) + v cilc 
H(z)r 



r < 10r v i r ; 
r ^ 10r v ir 



r vil < r < 10r vir 
r ^ lOrvir, 



(4) 



where H(z) is the Hubble constant at redshift z. 



3.2 Photoionisation of the IGM 

If the galaxy spectrum is of the form J(y) oc u 13 blueward of 
the hydrogen ionisation threshold, vh, then the total pho- 
toionisation rate at distance r from the galaxy is 



Lion = Bc(r)r BG + 



(3 - 3 47rr- 2 



(5) 



where we approximate the hydrogen photoionisation cross 
section as au(y) = &o(y /vh)~ 3 , in which oo = 6.3 x 10~ 18 
cm" 2 . The photoionisation rate due the extra-galactic back- 
ground, Fbg, has been determined observationally to be 
~ 10" 13 s" 1 between z = 5 and z = 6 (|Fan et al.| [2006). The 
factor B c (r) is the (radius dependent) boost in the ionising 
background due to the clustering of surrounding galaxies. 
Clustering is ignored in our fiducial model, i.e. B c (r) — 1, 
but may be as large as 10 2 (see S 13.3)1 . For a given Fi on , the 
neutral fraction of hydrogen is calculated from 



xh = 1 + 



Fio 



r r io 



/ r io 



2-1 1/2 



_|_ / x ion \ ' ^-g^ 
rcc ' J 

which was derived assuming photoionisation equilibrium and 
ignoring helium. Here, uh is the number density of hydrogen 
nuclei and the case-A recombination coefficient a rec = 4.2 x 
10 _13 (T gas /10 4 K) ' 7 cm 3 s" 1 (e.g. iHui fc Gnedinl Il997t ). 
We have verified that the total optical depth to ionising 
radiation is negligible in our calculations. We note that 
our optically-thin approximation is invalid if a dense, self- 
shielded and optically thick Lyman-limit system (LLs) lies 
along a given sightline. If anything, the presence of LLs 
would lower the actual photoionisation rate, which would 
increase the neutral fraction of hydrogen, which would in 
turn lower the overall Lya transmission. 

3.3 The Impact of Clustering of Sources on the 
Ionising Background 

We have calculated the enhancement of the ionising back- 
ground due to clustering of nearby galaxi es using a semi- 
analyt ic model based on that described in IWvithe fc Loebl 
(2005). Details of the model can be found in Appendix A, 
and we summarise only the its main features here. We define 
a mean free path of ionising photons by A, which we assume 
to be independent of radius. At a distance r from the central 
galaxy, this flux may be written as 



^b g (r) 



d 3 x- 



^(1*1) 



4tt| 



exp 



A 



(7) 



where dN/dt(r) gives the ionising photon production rate 
per unit volume from sources surrounding the central galaxy 
at radius r. The model computes dN/dt(r) using a radius 
dependent mass function for galaxies in the infall-region. 
More details are given in Appendix A (Eg IA1|) . Note that 
several model assumptions enter the expression of dN/dt(r), 
which introduces some uncertainty into the exact enhance- 
ment of the background flux. However, the main purpose of 
this calculation, is to show that the boost in the ionising 
background is large (this statement is independent of the 
detailed model assumptions). 
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10 ^ A-5 Mpc 





Radius (physical Mpc) 

Figure 2. The boost in the ionising background B c {r) due to 
clustered sources surrounding the central galaxy for the cases 
Aftot = lO 1:L Af0 {top panel) and Af to t = 1O 1O M {bottom panel). 
The boost is shown for 3 values of the assumed mean free path of 
the ionising photons (text). The vertical axis on the right shows 
the associated physical photoionisation rate (in s — 1 ). For com- 
parison, the photoionisation rate due to the central galaxy is 
shown as the thick grey line for Af* = 10 2 AfQ/yr {top panel) 
and Af, = 10AfQ/yr {bottom panel also in units of the measured 
Fbg)- For both cases, B c {r) -> 10 2 when A ^2 Mpc. Clust ermg 
of sources may boost the ionising background significantly. As 
a result, the total photoionisation rate is not dominated by the 
central galaxy beyond a few virial radii. 
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Figure 3. The impact of a post-reionisation i.e. fully ionised IGM 
on the observed shape of Lya emission line emitted by a z = 5.7 
star forming galaxy (Af* = lOAfQ/yr) in a dark matter halo with 
a mass of Aftot = 10 Mq (this mass is relevant for the extent 
of the infall region). The intrinsic shape {dashed line) of the Lya 
line is Gaussian. The solid represents the observed shape after 
IGM processing. Here, only ~ 12% of the Lya is transmitted to 
the observer. 



The resulting ionising background is shown in Fig- 
ure [2] for the cases in which the central galaxy sits in 
dark matter halos of mass M to t = 10 M@ (top panel) and 
Aftot = 1O 1O M {bottom panel). The flux F hg {r) is written 
in units of the mean background flux [i.e. we have plotted 
B c {r), Eq fl5J], To show the dependence of the boost on 
the assumed A, B c {r) is shown for A = 1,2 an d 5 physi- 
cal M pc, which covers the range of A derived bv lFan et al.l 
(2006). Note that B c (r) -» 10 2 for A ^ 2 Mpc in both cases. 
For comparison, the ionising flux from the central galaxy is 
shown as the thick grey line for M, = 10 2 Mq /yr {top panel) 
and M, = lOMe/yr {bottom panel, more precisely, F ga i/rBG 
is shown, see Eq[5]). In all cases, the total ionising flux is not 
dominated by that of the central galaxy beyond a few virial 
radii. The impact of this enhancement on the Lya line is 
discussed in § 14.11 

3.4 The IGM Transmission 

The total opacity seen by a photon initially at frequency v 
is 

t(v) = a f ds nn(s) xn(s) <j>(x(v) - mGM ^ \ (g\ 

Jr vil V Vth > 

where we denoted frequency in units of the dimensionless 
variable x = {v — v a )/Ai>D, in which vd = fa{vth/c); 
Vth = (2fcT g as/w,p) 1 ^ 2 is the thermal speed of atoms in 
the gas at temperature T ga s- We assume the gas temper- 
ature to be Tg as = 10 4 K, which is appropriate for the pho- 
toionised IGM. Note that the temperature of the IGM may 
actually have been as high as T gas = 5 x 10 4 K shortly a f- 
ter reionisation was completed (e.g. iHui fc Haimanl f2003). 
However, we found that our results are practically insen- 
sitive to the precise value of T gas . The Lya absorption 



cross section is written as o~(x) = o~o(j>(x), where o~a = 
5.67 x 10" 14 (T gas /10 4 K)~ 1/2 cm 2 and <j>{x) is given by 

*{*) = { ^^a COTe; ■ 0) 

1 ~ vk? wm s> 

The gas density profile nn{r) given by Eq @, is azimuthally 
averaged. In reality, for a given line of sight, the density 
fluctuates around this average. These density fluctuations 
can cause significant deviations from the opacity calculated 
from Eq (JHJ . If the density fluctuations are characterised by 
a distribution P{A)dA, which gives the probability that the 
gas density at radius r is in the range p(r)[A±dA/2] (where 
p{r) is given by Eq [4j| , then the average IGM transmission 
at frequency x 

/•A max 

(e" r )(x) = / dA P(A) e _A t(x) . (10) 
Jo 

This equation is equivalent to Eq (9) in Barkana & Loeb 
(2004). Here, we take A max = 100, but note that the 
results presented depend very weakly on its value. The 
function P(A)dA is not well known. Following Barkana & 
Loeb (2004) we take the d istribution P(A)dA derived by 
iMiralda-Escude et all (|2000T ) , who obtained this distribution 
at various redshifts from a simulation. This means that the 
density fluctuations in the overdense regions are assumed to 
be the same as those at mean density. 

Since 4>{x) is peaked strongly around x = 0, Eq (JSj 
states that the largest contribution to t{x) comes from gas 
that resonates with the photons at this frequency. By com- 
bining Eq ((51 , Eq © and Eq (J5J we can calculate the impact 
of the IGM on the Lya emission line. 
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4 RESULTS 

Our fiducial model consists of a star forming galaxy embed- 
ded in a dark matter halo with a total mass of Mtot = 
10 11 Mq. At z = 5.7, the halos' associated virial radius 
and circular velocity are 23 kpc and 133 km s _1 , re- 
spectively. The value of Mtot is chosen to facilitate com- 
parison with Santos (2004). However, this choice of Mtot 
may be justified by the observed bias and clustering of 
high-redshift LAEs which suggest M to t = 1 11 - 1O 13 M 
|Ouchill2003l y Furthermore. I Wyithe fc Loebl (|2006n have re- 
cently presented a simple model which reproduces t he ob - 
served Lya luminosity function of Kashika wa et alJ |2006). 
In their model, the observed LAEs have Mtot = 10 10 — 
lO n M 0. Using detaile d j>emi-analytic galaxy formation 
models ILe Delliou et al.l (|2006h found the typical halo mass 
of the observed LAEs to be M to t ~ lO n M . For the impor- 
tant results presented in this paper, our analysis is repeated 
for the case M to t = 1O 1O M0. For larger M to t, the infall re- 
gion is more extended and the IGM transmission is expected 
to be (even) lower. We assumed a star formation rate of 
M* = IOMq /yr, and the escape fraction of ionising photons 
to be /esc = 0.1. The spectrum of the galaxy is of the form 
J(v) oc u 13 , where /3 = —4.5, for v ^ uh (§ 12 . 3|) . The impact 
of clustered sources on the ionising background is discussed 
later [here, B c (r) = 1, Eq[5]. According to Eq (|H2[I . this 
galaxy has an intrinsic Lya luminosity of LL ya = 2 x 10 43 
erg s _1 . 

In Figure [3] we show the shape of the Lya line after 
it has been processed by the IGM. The solid/ dashed lines 
show the observed/intrinsic Lya? spectrum, respectively. As 
discussed in §[TJ the blue side of the line (x > 0) is suppressed 
significantly by the IGM. However, the suppression extends 
into the red part of the line. This is because of resonant ab- 
sorption by residual neutral hydrogen gas that is falling onto 
the galaxy. The peak of the observed spectrum is redshifted 
relative to the true line center by an amount set by the infall 
velocity. Here, the redshift is viGM(rVir) = — ^ch-c- 

We denote the flux density in the Lya line by J(x), 
where x is the dimensionless frequency variable introduced 
in § 13.41 The total transmission is then given by 

where (e~ T )(x) is given by Eq (|10|l . For the example shown 
in Figure[3] T a = 0.12, i.e. only 12% of the flux that emerges 
from the galaxy is transmitted through the IGM. The total 
observed flux of this galaxy is then T a xL-L ya /[4Tvd L (z)] — 
6 x 10~ 18 erg s" 1 cm" 2 . 

The exact shape of the observed Lya line and the total 
transmitted flux depends on several other parameters, such 
as the assumed width of the intrinsic width of the Lya line; 
the total ionising luminosity of the galaxy; the value of the 
ionising background etc. The variation in observable Lya 
properties with these, and other parameters is studied next. 

4.1 Variation of the Line Shape 

In Figure [4] we show the impact of various parameters on the 
shape of the Lya line after processing through the IGM. The 
IGM transmission for each curve is show n as a label. Also 
shown is the 'skewness parameter', Sw l|Kashikawa et all 



2006)0- This parameter quantifies the asymmetry of an emis- 
sion line, which can be used to separate high redshift LAEs 
from lower redshift interlopers. This separation is possible 
as absorption in the IGM causes high redshift Lya emission 
lines typically to be asymmetric (with a positive/negative 
skewness when the spectrum is shown as a function of wave- 
length/frequency, respectively), while lower redsh i ft Ha , 
[0111] or [OH] emitter s are not. Kashika wa et al.l (2006); 
IShimasaku et al.l l|2006h found that the requirement Sw ^ 3 
A can effectively distinguish between the two populations. 
Calculating Sw f° r our models is not trivial. The largest 
contribution to the skewness comes from pixels that lie far 
from the observed Lya peak flux, and we found that the 
precise value of Sw depends somewhat on the assumed flux 
threshold, /mm, below which the Lya line is lost in the noise. 
In particular, a very low level of residual flux on the blue 
side of the Lya line can cause Sw to be negative for some of 
our models. To account for the dependence of Sw on /min, 
we have shown the values of Sw we found for each model 
for the range 0.01/ max sS /mm =S 0.1/ max , where / max is the 
observed peak Lya flux density. 

• Top Left: Variation of M*: Here, we have plotted the 
cases M, = 10M Q /yr (black), M* = 10 2 M Q /yr (grey) and 
M* = lO 3 M0/yr (red). Increasing M* results in a higher 
output of ionising photons which reduces the neutral frac- 
tion in the infall region, and allows more flux, both red 
and blueward of the line, to be transmitted. Note that for 
M» > 100Af Q /yr) Sw < 3 A in difference to known high- 
redshift LAEs. For M* = lO 2 M0/yr the spectrum shows a 
long tail on the blue side of the line, which makes Sw nega- 
tive. On the other hand, for M* = lO 3 M0/yr the spectrum 
becomes practically symmetric with a skewness parameter 
-0.7 <S W <0.5 A. The precise value of M, at which Sw be- 
comes less than 3 A depends on the size of the infall region, 
and therefore on Mtot- For example, given M to t = 10 10 M Q , 
Sw <,3 A for M» ^10M e /yr. 

• Top Right: Variation of the width of the intrinsic Lya 
line: a a — 1.0u c irc (black), a a = 0.7u c irc (red) and a a — 
1.5v c irc (grey). Clearly, increasing o a increases T a , because 
a larger fraction of photons lie redward of the wavelength 
range affected by the infalling gas. Furthermore, the sharp 
break in the profile becomes more prominent towards larger 

Oct- 

• Lower Left: Variation of the ionising background Tbg: 
We show cases in which the ionising background is Tbg 
(black), IOFbg (grey) and IOOFbg (red). The main effect 
of raising the intensity of the ionising background is to in- 
crease the transmission on the blue side of the line. The 
skewness parameter, Sw, varies non-monotonically for sim- 
ilar reasons as those mentioned above: because of the long 
blue tail in the spectrum, Sw is very negative for the case 
of IOFbg , while the spectrum becomes almost symmetric 

2 Here, we denote the flux density in pixel number 'i' as f(xi). 
The parameter Sw> which has units of A, is defined as Sw = 
[A10 ,rcd -^10, blue] 5\ where AiQ re( j 

/blue are the observed wave- 
lengths at which the flux drops to 10% of its peak value, and S = 
J^J27=l( x -i -x) 3 f( x i)> where I = J2f( x i), x = J2 x if( x i)/ 1 
and a = [*%2(xi — x) 2 f(xi)/ 1] 1 ' 2 . For a detailed definition and 
discussion of the skewness parameter, t he reader is referred to 
the Appendix of Kashikaw a et all l l2006t) . 
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Figure 4. The observed Lya line from z = 5.7 galaxies in a reionised IGM for a range of different models. The upper right corner of 
each panel shows the model-parameter that is varied. Top Left: M* = IOMq/jt (black), Af» = 10 2 Mq /yr (grey) and M* = 10 3 MQ/yr 
(red). Top Right: cr a = 1.0u c j rc (black), a a = Q.7w c j rG (red) and cr a = 1 .5n c i rc (grey). Lower Left: Variation of the ionising background: no 
ionising background (dotted), VsQ(black), IOFbg (grey) and IOOFbg (red). Lower Right: Impact of galaxies peculiar velocity: i> pe c = 
(black), Vp C c = — 0.5<T a (red) and v pcc = 0.5o a (grey). For each model the total transmission, T a and the skewness of the line, Syy, are 
shown. 



again for the case of IOOFbg- The total IGM transmission 
depends only weakly on the intensity of the ionising back- 
ground. The reason is that the transmission is determined 
strongly by the neutral fraction in the infalling gas, which 
is set primarily by the strength of the ionising luminosity 
of the galaxy. This can be understood as follows: define 
the 'proximity radius', r prox , as the radius at which the 
photoionisation rate due to the galaxies ionising photons 
equals that of the ionising background (Fbg = TGaO- F° r 
the fiducial model, the proximity radius lies at r prox ~ 200 
kpc (~ 10r v i0, which is well outside the infall region. For 
r -C r pr0 x, the ionisation rate due to ionising background 
makes up only ~ (r/r plox ) 2 ~ [0.01 — 0.25] of the total ioni- 
sation rate within the infall region [in which r ~ (1 — > 5)r v i r , 
see § 13.1] . As a result, raising Fbg by a factor of ~ 10 only 
contributes significantly to F to t in the outer layers of the 
infall region. 

Also shown is the spectrum that would be observed when 
the ionising background is absent, Fbg = 0.0 (dotted line). 
The damping wing of the neutral IGM (outside of the HII 
sphere created by this galaxy), significantly reduces the flux 
redward of the infall feature as well. This damping wing 



optical depth is only impor tant when the neutral fraction in 
the IGM is large x m £0.1 (|Santosll200# l 

• Lower Right: Impact of galaxy peculiar velocity: v pcc = 
km s _1 (black), v pcc = — 0.5cr a (red) and v pcc = 0.5a a 
(grey). A negative peculiar velocity corresponds to the 
galaxy receding from the IGM. Clearly, more flux is trans- 
mitted towards larger negative velocities. The choice v pcc — 
—v c i IC is equivalent to models that use a static IGM, and 
would result in T a ~ 0.5. Furthermore, the shape of the 
Lya line is quite insensitive to v pac : all lines are asymmetric 
with a sharp cut-off on the blue side of the line. 

Lastly, we investigate is the impact of the radial dependent 
enhancement of the ionising background due to clustered 
sources. We showed in § !3.3l that this enhancement, denoted 
by B c (r), can be as large as 10 2 , particularly in the close 



3 The size of the Stromgren sphere was calculated from Qjj X 
tgal = J fls knr 2 nH (r)dr, in which t ga i = 3 X 10 s yr, is the life 
time of the galaxy. The neutral fraction in the IGM for r > R B 
as set to 1, while for r < R s photoionisation equilibrium was 
assumed (EqlgJ. 
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Figure 5. The impact of an enhanced ionising background due 
to clustered sources surrounding a galaxy on its observed Lya 
line. The red dashed line is the fiducial case with B c (r) = 1, 
whereas the black solid line shows the case in which B c (r) takes 
the form shown in Fig. [2] (the case of A = 1 Mpc). The blue 
side of the line is significantly enhanced, resulting in an increased 
transmission, T a = 0.30. Note that because of the skewness of 
— 13A $5 S w — 12A, this LAE would have been misidentified 
a low redshift interloper (if the Lya line were its only observed 
feature, see § 15.41 1. 



vicinity of the galaxy. In Figure [5] we show the observed 
Lya line for the profile B c (r) presented in the top panel of 
Fig. [2] (black solid line, X = 1 Mpc). The red dashed line 
represents the fiducial case with B c (r) = 1 for comparison. 
The figure shows that the blue side of the line is significantly 
enhanced (as we noticed above, see the lower left panel of 
Fig HI i resulting in a more prominent 'dip' in the spectrum. 
The overall transmission is enhanced by almost a factor of 
3 to T a — 0.30. The skewness of the line is — 13A ^ S w ^ 
— 12A, and therefore, this line would not be identified as 
a high- redshift LAE (see § 15.41 for further discussion). It 
is interesting to note that since A decreases with redshift, 
the boost in the ionising background increases with redshift, 
which in turn implies that more Lya lines may have Sw <^A 
towards higher redshift. 

We have found the observed shape of the Lycv line to 
vary widely, depending on M,, Tbg, cr a and the systemic ve- 
locity of the line, « poc - The main results can be summarised 
as follows: 

(i) For large star formation rates and strong boosts 
of the ionising background by clustering of nearby sources, 
the Lya line loses its trademark asymmetry (the skewness 
parameter for these models is Sw <J3 A) . This may be- 
come more common towards hig her redshift. iTaniguchi et al.l 
|2005l ); lKashikawa et all (|2006h found 20-30% of their LAE- 
candidates to have symmetric lines. It could be that at least 
some of these are LAEs, as opposed to [O II] A3727 emitters. 
This is discussed in more detail in § 15.41 

(ii) For most models, however, the Lycv line is asym- 
metric with a sharp cut-off on the blue side, and skewness 
parameters of Sw J>3 A. In such cases, the IGM scatters 
more than half of the Lya flux out of our line of sight, with 
mean transmissions of typically T a = 0.1 — 0.3. Further- 



more, infall causes the peak of the Lya line to be redshifted 
relative to the true line center. 

(ill) The Full Width at Half Maximum (FWHM) of the 
line is <,u c irc for most models, apart from those with large 
M t , ^pec, ®oi or S c (r). Models with large Tl^f* and/or S c (r) 
have Sw <3 A, and so produce either broad symmetric lines 
or lines that are asymmetric with skewness parameters that 
are negative. 



4.2 The Mean IGM Lya Transmission Around 
Galaxies 

The previous discussion focused on the change in the Lya 
line shape and Lya transmission as a function of star- 
formation rate, ionising background, the width of the in- 
trinsic Lya and the peculiar velocity v pcc . These calcula- 
tions focused on z = 5.7 galaxies. Here, we investigate the 
total transmitted flux for the fiducial model in an ionised 
IGM (i.e. /osc = 0.1, Vp CC = 0, a a = u circ and Z = 0.05) 
at additional redshifts of z = 4.5 and z = 6.5 as a function 
of M». Thus we have assumed that reionisation was com- 
pleted by z = 6.5, although evidenc e exists from high red- 
shift quasars may suggest otherw ise (|Wvithe fc Loebil2004bl ; 
iMesinger fc Ha iman 2004, 200(1). However, the analysis pre- 
sented here also applies when galaxies at z = 6.5 are embed- 
ded in large HII bubbles which are expected to exist prior 
to the end of reionisation (see § [S3] for a more detailed dis- 
cussion) . The impact of varying the other model parameters 
can be inferred from Figure [4] 

The results of this calculation are shown in the lower 
panel of Figure [U The black, red and grey lines represent 
galaxies at z = 4.5, z — 5.7 and z — 6.5, respectively. The 
grey dotted line represents the scenario in which the ionising 
background is set to at z = 6.5. The damping wing of 
hydrogen in the neutral IGM strongly suppresses the trans- 
mission as (we noticed already in § 14. 1[1 . At all redshifts, 
T a ~ 0.1 for M* <JlOM /yr. For larger star formation rates, 
the mean transmission increases slowly and reaches the value 
T a = 0.3 only when M» ~ 9OM /yr, M, ~ 3OOM /yr and 
M* ~ 5OOM /yr at z = 4.5, z = 5.7 and z = 6.5, re- 
spectively. Figure [6] strongly suggests that not properly ac- 
counting for scattering in the IGM, can lead to significant 
underestimates of Li. ya , and thus of the star formation rates. 
In the upper panel, we show the observable Lya flux as a 
function of star formation rate. The solid/dashed lines rep- 
resent calculations in which the transmission was/was-not 
accounted for. The thin horizontal lines denote the detec- 
tion thresholds of rec ently perfo rmed Lya surveys at z = 4.5 
(iDawson et al.ll2004h. z = 5.7 (|Shimasaku et al.l 120061 ) and 



6.5 |Taniguchi et alJliooBT) . As an example, for a given 



Lya flux from a 



5.7 galaxy, say S a = 10 erg s 



formation rate of M* = 2Mq /yr, while the actual star for- 
mation rate was M* ~ 2OM /yr. Therefore, not accounting 
for IGM scattering could lead to substantial underestimates 
of the Lya-derived star formation rates occurring in LAEs. 
We discuss this in more detail in S 15.11 
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Figure 6. Bottom Panel: Mean IGM transmission for a galaxy in a halo of mass Aftot = 1O 11 M0 as a function of the star formation rate 
(M*), at 2 = 4.5 (black), z = 5.7 (red) and z = 6.5(grey). The dotted grey line represents the scenario in which the IGM is completely 
neutral at z = 6.5. Top Panel: The observed Lya flux for our fiducial model as a function of M,, at the same redshifts presented in the 
bottom panel. The solid lines and dashed lines show the observed flux when the IGM transmission is/is not accounted for respectively. 
For all SFRs the mean transmission is significantly less than 1, which may lead to an underestimate of M*. The thin horizontal lines 
denote the detection thresholds of recently performed Lya surveys at z = 4.5, z = 5.7 and z = 6.5 (see text). 



4.3 Lya Transmission in the Absence of a Virial 
Shock 

Our model assumes that no absorption takes place inside 
the virial radius. The motivation for this assumption is that 
the virial temperature is well above 10 4 K for the halos con- 
sidered in this paper, and as a result, the virialised gas is 
collisionally ionised to a very high level (in our models the 
collisional ionisation rate of hydrogen in the virialised gas is 
~ 10 -10 s^ 1 ). However, the cooling time of this gas is short. 
In particular, the co oling time is shorter than the dynamical 
time of the system l|Haiman et al.ll2000h , and hence the gas 
rapidly loses its pressure support and collapses further. In- 
deed, it has been argued that because of this efficient cooling, 
the virial shock does not exist in halos with masses below 
M <,2 x lO n M (the precise mass depends on the metallic- 
ity of the gas, Birnboim & Dekel, 2003). We may therefore 
expect that infall continues at radii less than the virial ra- 
dius. The precise infall velocity profile of this gas is not well 
known. However, the absence of a virial shock would re- 
duce the Lya transmission further, especially in cases where 
the infall within the virial radius occurs at speeds exceed- 
ing t>circ, because then, the absorption would extend even 
further into the red part of the line. 



5 DISCUSSION 

5.1 Deriving Star formation Rates from Lya 
Fluxes 



formation rate. IShimasaku et all (|2006l ) found 28 spectro- 
scopically confirmed LAEs. They found that the Lya derived 
star-formation density for th eir survey volume lies a factor 
of ~ 3 below the estimates of lBouwens et all ll2006l, down to 
the same UV- luminosity). Recently, iTaniguchi et ail (|2005l ) 
found that the Lya derived star formation rates lie, on av- 
erage, a factor of ~ 5 below values derived from the UV- 
continuum measurements. The typical UV-derived star for- 
mation rates for galaxies in their sample were in the range 
~ 5 — 45Mq /yr, which according to our Figure [6] would 
translate to T a ~ 0.1 — 0.2. Scattering in the IGM therefore 
fully accounts for the apparent discrepancy between UV and 
Lya derived star formation rates. 



In § 14. 21 we have shown how Lya scattering in the IGM can 
reduce the detected Lya flux by up to an order of magni- 
tude. This implies that estimates of the star formation rate 
based on Lya fluxes generally underestimate the true star 



Note however, that the discrepancy between UV and 
Lya derived star formation rates, is not considered a seri- 
ous issue. The reason for this is that Lya photons may be 
subject to a different amount of dust attenuation than UV- 
continuum photons. For example, scattering of Lya photons 
through the ISM increases their total traversed path through 
the galaxy, making Lya photons more vulnerable to dust 
attenuation (on the other hand, in the case of a multiphase 
ISM, the Lya may be subject to less attenuation, see § 12. 3p . 
Therefore, the observed discrepancy between Lya and UV- 
derived star formation rates can also be interpreted as result- 
ing from Lya b eing preferentially destroyed by dust inside 
the galaxy (see iKunth et al1 IT998. for a discussion on the 
difficulties associated with estimating star formation rates 
from Lya fluxes). However, it is interesting that our low 
value of the transmission through an ionised IGM provides 
a natural explanation for the discrepancy. 
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Constraints from the Observed Lya Line 
Width+EW. 



5.2 



Figures [3151 show that the typical observed line width of our 
model LAEs is ~ 2 — 10 A, with the fiducial model hav- 
ing a FWHM of ~ 3 A . The bulk of the observed LAEs in 
IShimasaku et alJ l|2006l ) have a FWHM in the range ~ 2-14 
A, with the majority of the LAEs having FHWM ^5A. 
These values of the FWHM can be explained with larger 
v pcc or (J a . A larger star formation rate, or a boost in the 
ionising background have the same effect, but probably not 
while maintaining S w :>3 A. In particular, the high end of 
the range of observed FWHM cannot be explained this way. 
An elegant explanation for the large line width systems is 
provided by back-scattering off a superw ind driven shell of 
neutral gas l|Ahn et al. 2003; lAhn|[200l ). Here, Lya emit- 
ted by the galaxy scatters off an thick neutral shell that 
is receding from the galaxy with a velocity i> cxp . Photons 
that are scattered back to the galaxy have now picked up a 
redshift of 2« oxp , an d may make it th rough the shell with- 
out being scattered. lAhn et all (|2003l ) have shown that the 
Lya line can easily be broadened by several times u cxp . This 
back-scattered Lya would all be on the red-side of the line 
cent er, and would thus be t ransmitted through the IGM. 

iKashikawa et all l|2006h have found that marginal evi- 
dence exists for an anti-correlation between the Lya FWHM 
and the derived Lya luminosity, i.e. the measured Lya flux. 
This is contrary to what one would expect for the superwind 
scenario sketched above (in which the IGM transmission 
is larger) unless the wind-driven shells contain dust which 
could reduce the escape fraction of Lya photons from the 
galaxy. An alternative explanation for this observation may 
be that that the broadest lines are broadened by resonant 
scattering (although we find this more unlikely as discussed 
in 6B. 



5.3 Comparing the Lya and LGB Luminosity 
Functions. 



IKashikawa et all ()2006h compared the restframe-UV lumi- 
nosity function of their z = 6.5 LAE sample w ith that of 
the i-drop-out sample of iBouwens et all (120061) and found 
that the faint end slop e for LAEs at z = 5.7 jHu et al.|[2004l ; 
IShimasaku et al.ll2006h and z = 6.5 l|Kashikawa et alj|2006h 
is flatter than that of the i-drop out galaxies. The differences 
between the LAE and i-drop out UV-lu minosity functions 
only b ecomes prominent for Muv > — 20.5. IKashikawa et ail 
l|2006t ) suggest that this is most likely due to a combina- 
tion of detection incompleteness at faint Lya and UV-fluxes 
(Muv Jj — 20.3 corresponds to the 3 — a detection threshold). 
Using the standard conversion from UV-conti nuum flux den- 
sity to star formation rate fe.g. lKennicuttll998l ). we find that 
M uv ^ - 20.5 corresponds to M, <; 8.OM yr" 1 . Accord- 
ing to Figure [6] the IGM transmission associated with these 
star formation rates is low, T a <^0.12 in this case, and the 
expected Lya flux is S a <J6 x 10 -18 erg s _1 cm -2 ; almost ex- 
actly the detection threshold. For lower star formation rates 
(and higher Muv) a detection in Lya requires a larger IGM 
transmission around these galaxies, which requires a large 
Wpec, or a a (§ 14.111 . This suggests that at fainter Lya (and 
UV) fluxes indeed only a small subset of LAEs is currently 
detected. 



5.4 Symmetric and Double Peaked Lya Emission 
Lines; A Source of Selection Bias? 

In § !4.1l we found that the Lya line becomes either symmet- 
ric or asymmetric in the wrong direction (Sw <J3 A) for large 
star formation rates and for strong enhancements of the ion- 
ising background. This suggests that the current criterion of 
asymmetry introduces biases against detection of galaxies 
undergoing intense episodes of star formation and/or em- 
bedded in regions of enhanced ionising background. In prac- 
tice however, the criterion Sw > 3A is used in combina- 
tion with the broad band p hotometry in the B and V-bands 
(e.g. IShimasaku et alj|200q ). High- redshift LAEs are practi- 
cally invisible in both these bands due to absorption in the 
IGM blueward of the rest-frame Lya frequency. Therefore, 
only those objects having both Sw < 3A and a detection 
in the B or V Band are classified as low redshift interlop- 
ers. Objects for which Sw < 3A that are not detected in B 
and V are classified as an "unc lear" objects. Approxim ately 
15% of the objec t s dete cted by IShimasaku et all |2006f) and 
IKashikawa et al.l {2006) were classified as such and could 
thus be high -redshift LAEs. S i x out of 34 of the S w > 3 A 
objects from IShimasaku et alJ (|2006h are bright in B or V. 
IShimasaku et al.l ~i 20061 ) argue that these objects are gen- 
uine LAEs, with the B/V flux contributed by nearby faint, 
contaminating sources. Likewise, some of the 19 objects in 
their 63-object spectroscopic sample with symmetric lines 
and detection in B/V could be LAEs and foreground B/V 
sources. However, the above considerations suggest a small 
fraction of Sw < 3A LAEs and argue against a serious se- 
lection bias. 

The low number of high-redshift LAEs with Sw <>3 A 
could be attributed to a small number of LAEs being embed- 
ded in regions of strongly enhanced ionising backgrounds, or 
to a small number of galaxies undergoing intense episodes 
of star formation. It could also imply that the IGM sur- 
rounding LAEs is denser than we have assumed. To make 
observed Lya lines have S w > 3 A, one needs to enhance 
the absorption in the IGM. Conversely, if we ignore infall, 
or lower the IGM density, then lower star formation rates 
or smaller enhancements of the ionising background are re- 
quired to produce Sw <J3 A-lines. We point out that al- 
though high-redshift LAEs with symmetric lines appear to 
be rare, these objects may be cosmolo gically very int erest- 
ing. For example, recent simulations bv lLi et all (|2006l ) show 
that the star formation rates in progenitor galaxies of the 
luminous z — 6 quasars may be as large as M* = 1O 4 M0 
yr -1 (at z — 6 — 14). Absorption in the IGM would almost 
certainly not produce Sw >3 A Lya lines. 

Some Lya emission lines shown in Figures HJ] and [5] 
with Sw <J3 A have two peaks that are separated by ~ 7 
A. Coincidentally, this is almost exactly the separation ex- 
pected for a resolved [Oil] doublet (~ 6.6 A , in the ob- 
served frame at z = 1.2). In cases where neither B nor V 
is detected, such a Lya emitter may be confused with a 
2 = 1.2 [O II] doublet. Whether or not a Lya line con- 
sists of two peaks depends on the following considerations. 
The largest contribution to the Lya opacity at a frequency 
x comes from the resonant region (defined as the region 
where the photons are at resonance, Eq[SJ). Therefore, if the 
resonant region lies a distance r from the central galaxies, 
the total opacity to photons at frequency x can be approx- 
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imated by t(x) ~ nn (r)xH(r)v t h/[dv/dr] oc n 2 H (r) / J(r). 
Here, we used the relation Xh(t) oc nu(r) / J(r), and as- 
sumed for simplicity that the variation of dv/dr with radius 
can be ignored. In our model, photons just blueward of the 
observed (red) peak, resonate with gas near the virial radius 
that falls inward at the maximum speed. If nii(r) decreases 
faster than [J(r)] 1 ' 2 with increasing r, then the opacity de- 
creases towards bluer wavelengths and produces a 'dip' in 
the observed spectrum. In our models nn(r) oc r _1 , while 
to produce a dip requires J(r) to change more slowly than 
r~ 2 , which is indeed the case when the ionising background 
is enhanced (Fig The precise location of the dip, and 
thus the separation of the two peaks, depends on the exact 
radial dependence of the density and velocity profiles of the 
surrounding gas, in addition to the ionising radiation field. 
Note that double peake d, spatially extended , Lya emission 
has b een observed (e.g. Ivan Oiik et al.lll997l : IWilman et al.l 
2005), often around radio galaxies. Various other mecha- 
nisms including superwinds and resonant s cattering can cre- 
ate d ouble peaked Lya emission lines (see IVerhamme et al.l 
2006, for a more detailed discussion on this issue). 



5.5 Lya Transmission at Higher Redshifts: z ^ 7 

Our analysis has focused on the redshifts bins z — 4.5, 5.7 
and z — 6.5, corresponding to windows of reduced atmo- 
spheric OH emission. Other win dows exist at z ~ 7, z = 8 
and z = 8.8 (Nilsson et al. 200<1). Lya emiss ion from higher 
redshifts is expected to be detected soon ( see live et al . 2006, 
for a detection of a z — 7.0 galaxy). Our analysis is ex- 
pected to be applicable here as well, even if reionisation was 
not completed until z — 6 — 6.5. The reason is that prior to 
overlap, large bubbles (with typical sizes of ~ 1 — 10 physical 
Mpc, hereafter pMpc) of ionised gas are expected to popu- 
late the I GM, each of which may h arbor many Lya emitti ng 
galaxies l|Furlanetto et al.ll2004alE; IWvithe fc Loebll2005h . 

The analysis presented by Santos! ( 20041 ) in which a 
galaxy sits inside its own HII region embedded in a neu- 
tral IGM, does not apply to most galaxies in this scenario. 
In particular, absorption by the damping wing of neutral 
intergalactic hydrogen is less important for galaxies inside 
these large HII regions. This is illustrated in Figure [7] Here, 
we show the Lya emission line of a galaxy at z — 8 (the other 
model parameters, except for Tbg, correspond to the fidu- 
cial model which is described in §2} for the case in which the 
galaxy is surrounded by its own HII region (which is ~ 0.7 
pMpc for our model), and for cases in which the galaxy is 
embedded within a larger HII region created by an ensemble 
of nearby sources. The IGM outside the bubble is assumed 
to be neutral. The other cases shown are for J?hii = 2.0 
pMpc (red line) and Rhu = 10.0 pMpc (grey Zine)Q 

Figure [7] shows that damping wing absorption on the 
red side of the Lya line center becomes decreasingly impor- 
tant towards larger Run- However, resonant absorption the 



x=(v-v n )/Av r , 
20 0° D 



-20 



4 We have assumed that the IGM is completely neutral outside 
the HII bub bles. Part i al ionisation of th e surrounding IGM (e.g. 
by X-Ravs. lohl l200ll; [Venkatesan et alj|200ll) . and neighboring 
bubbles reduce the contribution of the damping wing to the IGM 
transmission, in which case the dependence of the IGM transmis- 
sion on HII bubble size is even smaller. 
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Figure 7. The observed Lya line of a z = 8.0 galaxy (fiducial 
model parameters) for the case in which the galaxy carves its 
own HII region out of a fully neutral IGM (black line); when this 
galaxy lies in the center of a larger HII bubble of 2.0 pMpc (red 
line); and 10.0 pMpc (grey line). For each model T a and Sw> are 
shown. The presence of large 'uber-HII region' boosts the Lya 
transmission on the red side of the line. The dashed grey line 
represents the intrinsic spectrum. 



in infall region erases the Lya line for A <^1.095 x 10 4 A 
in all cases. Since the transmission is regulated primarily 
by resonant absorption in the infall region and by damping 
wing absorption in the neutral IGM, T a is practically insen- 
sitive to the neutral fraction of hydrogen inside the bubble 
at r ;>3 — 5r v i r • The overall transmission for the three cases 
is T a = 0.03 (Rmi = 0.7 pMpc), T a = 0.06 (Run = 2.0 
pMpc) and T a = 0.10 (Rmi = 10.0 pMpc). Therefore, at 
first glance, the presence of large HII bubbles prior to overlap 
boosts the Lya transmission by only a factor of ~ 3. We ob- 
tained similar results for a halo mass ofMtot = 10 10 M o ,for 
cases in which the Lya line is intrinsically narrower/broader 
(a a = 0.7/1.5-Ucirc), and for cases in which the galaxy pos- 
sesses a peculiar velocity with respect to the surrounding 
IGM (w P 8c = ±0.5u c irc). However, we caution that when 
calculating the Stromgren radius, we (incorrectly) ignored 
recombinations in the HII region. Therefore, the Stromgren 
radius we calculate for our model is actually too large (al- 
though note that ionising radiation produced by a cluster 
of fainter, undetected, sources surrounding the galaxy may 
increase the Stromgren radius, e.g. Wyithe & Loeb, 2005). 
If we decrease the Stromgren radius by a factor of ~ 3, then 
T a decreases by an additional factor of almost 2. Further- 
more, for galaxies in large bubbles the transmission could be 
boosted significantly due to nearby clustered sources (but 
note that this would likely result in Lya lines with Sw < 3 
A, see §E3J. 

The average HH-bubble size is related to the globally 
averaged ionised fracti on of the universe by volume, x- u y 
jFurlanetto et alj|200d ). and therefore it is useful to calcu- 
late the IGM transmission as a function of HH-bubble size 
for various models. The result of this calculation is shown 
in Figure [8] where we considered cases at z = 6.5 (dotted 
lines) and z = 8.0 (solid lines) for M» = IOMq yr" 1 and 
M t — 1O 2 M0 yr _1 . The minimum radius for each line cor- 
responds to cases of isolated galaxies in an otherwise neu- 
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Figure 8. The fraction of Lya that is transmitted through the 
IGM as a function of Hll-bubble size prior to the completion of 
reionisation. We have shown cases at z = 6.5 (dotted lines) and 
2 = 8.0 (solid lines) for M„ = 10M Q /yr and M, = lO 2 Af /yr. 
The minimum radius for each line corresponds to cases of iso- 
lated galaxies in an otherwise neutral IGM. The globally aver- 
aged ionised fraction of the universe by volume, ffj.y, is re lated 
to the characteristic Hll-bubble size (Furlanctto ct al. 2006). An 
approximate relation between i?Hll an d %i,v ls a l so shown. 



tral IGM. The relation between HII bubble size and neu- 
tral fraction depends on the nature of the ionising sources, 
the abundance of Lyman-limit systems, and o n the redshift 
|Furlanetto et al.ll200rj ; iMcQuinn et al.l l2006t ). An approx- 
imate relation between t he characteristic HII b ubble size 
Rmi and Xi t y derived by iFurlanetto et al.l (|2006h is shown 
in the Figure for 0.4 < xi )V < 0.9. The figure shows that the 
transmission drops only by a factor of ~ 3 compared to cases 
in which the IGM is fully reionised. As was mentioned above, 
part of the reason for this low contrast is that by ignoring 
recombinations within the Stromgren-sphere, we may have 
overestimated its size. Furthermore, clustered sources inside 
the large uber-HII bubbles could significantly boost the local 
ionising background, which would increase the transmission 
(see §E3J. 

This result has the following implication: as reionisation 
of the IGM proceeds, HII bubbles grow larger in time un- 
til they overlap. The transmission of Lya photons emitted 
by a galaxy inside such an HII bubble c ontinuously grows 
with bubble size l|Furlanetto et all l2006t ) . The precise de- 
pendence of T a on HII bubble size depends on the strength 
of the ionising background inside the bubbles . At overlap a 
'jump' in the ionising background may occur l|Gnedinll2000l. 
although see Furlanetto & Oh, 2005). If such a jump occurs 
it will not result in a jump in the IGM transmission around 
LAEs, as we showed in § 14.11 Therefore, LAEs of compara- 
ble luminosity should not appear to be significantly dimmer 
prior to overlap. The evolution in the Lya luminosity func- 
tion b etween z — 5.7 and z — 6.5 found bv lKashikawa et al.l 
(2006) is therefore not necessarily due to reionisation (unless 
the characteristic bubble size evolved very rapidly between 
z = 5.7 and z — 6.5, see Dijkstra et al. 2007 for a more 
detailed discussion). 



6 CONCLUSIONS 

We have calculated the impact of an ionised intergalactic 
medium on the Lya line emitted by galaxies at z 4.5. Our 
model accounts for gas clumping in the IGM and for the 
fact that high-redshift galaxies typically reside in overdense 
regions, which causes the velocity field of the IGM to deviate 
from the Hubble flow. We calculated the observed shape of 
the Lya line as a function of the galaxies ionising luminosity, 
the intrinsic width and systemic velocity of the Lya line and 
the extra-galactic UV-background (which may be enhanced 
due to unseen clustered sources by a factor of ~ 100). 

We found that the observed shape of the Lya line can 
display a wide range of shapes, depending on the precise 
values of A/*, Tbg, ""a and v pcc (see § 14.11 and Fig3|. Our 
most important findings are: 

• For large star formation rates, the strong asymmetry 
of the Lya line (believed to be the defining characteristic 
of a high-redshift Lya line) vanishes. For example, given 
Mtot = W^Mq (Mq = 1O 1O M ), a star formation rate of 
M, = 100M Q yr" 1 (M* = 10M Q yr" 1 ) produces S w < 0, 
whereas the value Sw ^ 3 A has been shown to distin- 
guish between confirmed high-redshift LAEs and lower red- 
shift interlopers. Also, we find that S w <>3 A for galax- 
ies in which the ionising background was boosted signifi- 
cantly (which may become more common towards higher 
redshift) by surrounding clustered galaxies. An enhanced 
ionising background may occasionally result in a double- 
peaked Lya emission line, mimicking a lower redshift [Oil] 
emission line. The criterion Sw 3 A could therefore intro- 
duce a selection bias against galaxies with high star forma- 
tion rates and/or galaxies embedded in an enhanced extra- 
galactic UV-background. However, at z=5.7 and z=6.5 this 
is not a serious issue as only ~ 15% of the known Lya emitter 
candidates has Sw <J3 A. Nevertheless, recent cosmological 
simulations have shown that star formation rates as large a s 
M* = 1O 4 M yr -1 may occur at z = 6- 14 (|Li et alj|2006h . 
Absorption in the IGM around these galaxies, while rare, 
would almost certainly not produce asymmetric Lya lines. 

• In most models, the Lya line is asymmetric (S w > 3 A) 
with a sharp cut-off on the blue side. In these models, the 
IGM typically scatters more than half of the Lya flux out of 
the line of sight, with a mean transmission of T a — 0.1 — 0.3. 

• The low values of the IGM transmission found in this 
paper suggest that not accounting for IGM scattering could 
lead to substantial underestimates of the Lya-derived star 
formation rates in LAEs. Although both the UV and Lya- 
derived star formation rates are subject to large uncertain- 
ties, our model naturally explains the apparent discrepancy 
that currently exists in the literature between the two. 

• Clustering of sources prior to overlap boosts the de- 
tectability of LAEs, as ensembles of clustered sources pro- 
duce large HII bubbles. This r educes the impact the damp- 
ing wing of the neutral IGM (jFurlanetto et al J 12006). The 
precise dependence of T a on HII bubble size depends on 
the strength of the ionising background inside the bubbles. 
Furthermore, the possible boost in the ionising background 
that follows the overlap of HII bubbles (when reionisation is 
completed) barely increases the IGM transmission. Increas- 
ing the local ionising background by a factor of ~ 10 only 
increases the IGM transmission by a few per cent. This sug- 
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gests that LAEs of comparable luminosity should not appear 
to be significantly dimmer prior to overlap. 

• The observed physical width of the Lya line contains 
useful information. In particular, Our models have difficul- 
ties reproducing observed FWHM >10 A , while conserving 
the asymmetry of the line, and we conclude that these LAEs 
are most likely to be galaxies with a superwind-driven out- 
flow. 

Our result that LAEs are not expected to become sig- 
nificantly brighter in a post-overlap IGM, naively implies 
that LAEs do not provide a direct probe of the end of reion- 
isation. However, our work suggests that LAEs are not ex- 
pected to dim significantly faster than oc dj / 2 (z) at z Jj7. 
We therefore expect Lya surveys to be very fruitful in un- 
covering young galaxies at these higher redshifts, which will 
allow studies of star formation in the very first galaxies. Con- 
straints on reionisation can then be obtained through var- 
ious methods which include for example, cross-correlating 
the distribution of LA Es with intergalactic 21-cm emission 
i|Wvithe fc Loebl I 20061 ); or by adding up the ionised vol- 
ume s urrounding individual LAEs as in Ma lhotra fc Rhoadsl 
(2006). Moreover, with a large enough sample of LAEs at 
various redshift bins, it is possible to detect a weak evo- 
lution of the IGM transmission, which is expected to ex- 
ist in our model due to the gradual growth of HII bubbles 
dFurlanetto et al.ll2006l ; IMcQuinn et al.ll2006l ; iDiikstra et~ai1 
12007 ) . The constraints on the IGM transmission can be fur- 
ther improved when Lya luminosity functions are combined 
with t he Lya emitters' rest-f rame UV luminosity functions 
(as in lKashikawa et al.ll2006l . also see Dijkstra et al. 2007). 
Additionally, since the abundance of observable emitters 
is modulated by the presence/absence of large HII bub- 
bles, the two-point correlation function of LAEs may be 
used to constrain the dependence of T a on H H-bubble size 
|Furlanetto et al.l l200rj ; IMcQuinn et al.ll2006l ), which puts 
constraints on the strength of the ionising background in- 
side the bubbles. 

Finally, we point out that our focus has been on the 
mean Lya profile. In reality, fluctuations in the transmis- 
sion from source to source are likely to exist. For example, 
the transmission is higher for sources viewed along a partic- 
ularly low density path. It is highly unlikely that the ma- 
jority of known LAEs are observed along such paths, as a 
high transmission is likely to be accompanied with Sw ^ 3 
A. In future work we will address this scatter in the Lya 
transmission in more detail using cosmological simulations. 
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APPENDIX A: THE BOOST IN THE IONISING 
BACKGROUND DUE TO CLUSTERING OF 
SOURCES. 

We model the ionising background in the overdense regions 
surrounding massive galaxies using a semi-analytic model for 
biased galaxy formation near high redshift galaxies. Since 
we focus on z > 4, beyond the peak redshift for quasar 
activity, we consider the case where AGN do not contribute 
to the back-ground. Making the assumption that the ionising 
radiation generated by galaxies was emitted over a Hubble 
time, and assuming the ionising photon mean-free-path (A) 
to be independent of density, we may derive an expression for 
the ionising background flux at radius R. We first calculate 
the ionising photon production rate per unit volume from 
sources surrounding the central galaxies 

Sra-jC"^)*^' ,AI) 

where nn is the hydrogen number density. In this expres- 
sion the number density of halos surrounding the central 
halo dn(R, Mh sio)/dM is found from the Press-Schechter 
mass f unction llPress fc Schechter! Ii"974h with the modifica- 
tion of lSheth fc Tormen ( 19991 ). expressed in units of (phys- 



ical Mpc) Mq . In the vicinity of a massive halo dn/dM 
must be modified relative to the background universe. The 
dominant contribution to the modification of the halo mass 
function originates from the accelerated formation of halos 
in overdense regions. T he calculation of this en hanced mass- 
function is described in IWvithe fc Loebl l|2005l ). 

The quantity D 4 in equation (|A1|I is the characteristic 
physical size of a spherical HII region generated by stars in 
a galaxy with a dark mat ter halo of mass A/h a i , and was 
given bv lLoeb et ail i|2005T ) 
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where /* is the fraction of baryons in the halo that are turned 
into stars, and / csc is the fraction of ionising photons that 
escape from the galaxy into the surrounding IGM. The num- 
ber of ionising phot ons per baryon inco rporated into Pop-II 
stars is A* ~ 4000 l|Bromm et alJ lioOl). and A rc i on denotes 
the number of photons required per baryon for the reionisa- 
tion of the HII region. The value of iV re i on depends on the 
number of recombinations and hence on the dumpiness of 
the IGM. 

When ( equation IA1[) is substituted into Eq (0, the en- 
hancement of ionising background at radius R can be com- 
puted. Several model parameters need to specified: First 
of all, we take Af m i n to correspond to w c i rc = 50 km s , 
which reflects the fact that only galaxies more massive than 
this can form in an ionised IGM, and hence contribute to 
the ionising background. Furthermore, we have assumed the 
ionising photons' mean-free-path to be equal to the mea- 
sured value at z — 5.7, A ~ 1 — 5 physical Mpc (|Fan et al.l 
2006). Last, we need to specify the precise dependence of 
D» on halo mass, by specifying the product /*/ eS c/A re i on . 
These model parameters were c hosen to achieve reionisation 
at z = 6 as in lWvithe fc Loebl (|2005l '>. 
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